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ABSTRACT 



Context. The putative tori surrounding the accretion disks of active galactic nuclei (AGNs) play a fundamental role in the unification 
scheme of AGNs. Infrared long-baseline interferometry allows us to study the inner dust distribution in AGNs with unprecedented 
spatial resolution over a wide infrared wavelength range. 

Aims. Near- and mid-infrared interferometry is used to investigate the milli-arcsecond-scale dust distribution in the type 1.5 Seyfert 
nucleus of NGC 3783. 

Methods. We observed NGC 3783 with the VLTI/AMBER instrument in the ^f-band and compared our observations with models. 
Results. From the A'-band observations, we derive a ring-fit torus radius of 0.74 ± 0.23 mas or 0.16 ± 0.05 pc. We compare this size 
with infrared interferometric observations of other AGNs and UV/optical-infrared reverberation measurements. For the interpretation 
of our observations, we simultaneously model our near- and mid-infrared visibilities and the SED with a temperature/density-gradient 
model including an additional inner hot 1400 K ring component. 
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The putative tori surrounding the accretion disks of ac- 
tive galactic nuclei (AGNs) play a fundamental role in 
the unification scheme of AGNs ( Anton ucci 1993). These 
tori most likely provide the material reservoir t hat feeds 
the accretion disk (e.g., ' Krolik & BegelmanI Il988h . Infrared 
interferometric observations can resolv e these mas-scale 
tori in the nea r-infrared (NIR) (e.g.. Swain et al. 2003; 
Wittkowski et al.ll2004: Kishimoto et al. : 2009a: Pott et al. 20ia 



Kishinioto et alj201 laE and mid-i n frared (MIR) (e .g., Jaffe et all 



2004"; 'Meisenheimer et al.' '20071; iTristram et al.' '2007', '2009', 
Beckert et al. 2008; Kishimoto et al. 2009b; Raban et al. 2009; 
Burtscher et all 12009. 2010; iTristram & Schartmannll201 11 and 



iKishimoto et alj 1201 lb . Paper I). In this paper, we present 
the first NIR interferometric observation of an AGN with the 
AMBERA^LTI instrument. 



2. Observations and data reduction 

We observed the Seyfert 1.5 AGN NGC 3783 in 2009 and 201 1 
(IDs 083.B-0212 an d 087.B-0578) with the ESO VLTI and the 
AMBER instrument dPetrov et al.ll200 7^. For these observations 
in the /T-band (see Table [1]), the AMBER low spectral reso- 
lution mode (LR) was employed. Figure [T] (top) presents ex- 



Send offprint requests to: G. Weigelt 
email: weigelt@mpifr.de 

* Based on observations made with ESO telescopes at Paranal 
Observatory under programme IDs 083.B-0212(A) and 087.B-0578(A). 



amples of target interferograms to illustrate the noise problem 
(K ~ 10.1). Long detector integration times (DIT) of 400 and 
800 ms were chosen to be able to recognize the faint fringes dur- 
ing data recording and correct drifts of the optical path dilfer- 
ences (OPDs) between the telescope beams. The interferograms 
shown in Fig. [1] are two-telescope interferograms, not AMBER- 
standard three-telescope interferograms. In 2009 (see Table [U, 
we recorded two-telescope interferograms since it is easier to 
correct OPD drifts in two-telescope than in three-telescope in- 
terferograms. In 2011, we recorded three-telescope interfero- 
grams, which provide closure phases. For data reduction of the 
two-telescope interferograms, we used our own software (de- 
veloped by one of us, KHH), which is able to reduce the non- 
standard two-teles cope interferograms. It is based on th e same 
P2VM algorithm (iTatulU et al.ll2007t IChelli et al.ll2009h as the 
AMBER amdlil^ software. To reduce the effect of the instanta- 
neous OPDs on the visibility, we applied a preprocessing method 
that equalizes t he OPD histograms of the target and calibrator 
interferograms (iKreplin et al. Figure [Tjand Table [T]show 

the derived visibilities, which are wavelength averages over the 
wavelength range of 2.0-2.40 jum. For data reduction of the 201 1 
three-telescope data, we used the standard AMBER data reduc- 
tion package amdlib version 3.0. Figure [1] (bottom) shows the 
closure phases of NGC 3783 derived from the three-telescope 
2011 interferograms. The average closure phase is 3.3 ± 26°. 
Closure phases are a measure of asymmetry. However, the large 



' The AMBER - reduction package amdlib is available at: 
http://www.jmmc.fr/data-processing_amberhtmi 
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Table 1. Observation log of our AMBER LR observations of NGC 3783 and its calibrator CD-37 739 1 . The data in the first 1 1 lines 
were observed in the two-telescope mode and the data in the last 4 lines in the three-telescope mode. The table lists the names, times 
of observations, projected baseline lengths, position angles PA, detector integration times DIT, seeing, number of interferograms, 
and derived target visibilities (errors +0.09). 



Name 


date 


Time of 


Telescopes/ 




PA 


DIT 


Seeing 


Number 


Target 






observation 


proj. baseline 




(°) 


(ms) 


(arcsec) 


of frames 


visibility 






(UTC) 


lengths (m) 














NGC 3783 


09/04/14 


01:19 - 01:31 


UT2-3/46.6 




28.8 


800 


0.93 


11x70 


0.89 


CD-37 7391 


09/04/14 


01:45 - 01:55 


UT2-3 






800 


0.84 


10x70 




NGC 3783 


09/04/14 


03:16-03:29 


UT2-3/45.2 




44.8 


800 


0.82 


11x70 


0.96 


NGC 3783 


09/04/14 


03:32-03:41 


UT2-3/44.8 




46.3 


400 


0.79 


10x70 


0.93 


CD-37 7391 


09/04/14 


04:31 -04:34 


UT2-3 






400 


0.78 


4x120 




NGC 3783 


09/04/14 


03:52 - 04:03 


UT3-4/62.0 




120.9 


800 


0.68 


10x70 


0.97 


CD-37 7391 


09/04/14 


02:11 -02:16 


UT3-4 






800 


1.45 


5x70 




CD-37 7391 


09/04/14 


02:00 - 02:05 


UT2-4 






800 


0.99 


5x120 




NGC 3783 


09/04/14 


03:44 - 03:49 


UT2-4/87.1 




90.0 


800 


0.72 


5x70 


0.86 


NGC 3783 


09/04/14 


04:12-04:20 


UT2-4/84.6 




94.7 


800 


0.81 


8x70 


0.85 


CD-37 7391 


09/04/14 


04:41 - 04:45 


UT2-4 






800 


0.87 


3x120 




NGC 3783 


11/05/18 


02:39 - 02:44 


UTl-2-4/51.2/113.8/ 


S0.3 


38.8/77.3/100.7 


400 


0.63 


7x120 




NGC 3783 


11/05/18 


02:50 - 02:55 


UTl-2-4/50.6/1 11.2/78.7 


39.8/78.9/102.8 


400 


0.68 


7x120 




CD-37 7391 


11/05/18 


02:00 - 02:06 


UTl-2-4/ 






400 


0.75 


7x120 




CD-37 7391 


11/05/18 


02:09-02:17 


UT 1-2-4/ 






400 


0.63 


7x120 





3. Geometric model fits 



Torus radius = 0.74 mas (21% AD flux, 0.5% host galaxy) 
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Fig. 1. Top: Examples of LR A^-band AMBER interferograms of 
NGC 3783 (~2.0-2.4/jm from bottom to top), which illustrate the 
noise problem. From left to right, the first and second interferogram 
are recorded with UT2-UT3 (46.6 m projected baseline) with DIT = 
800 ms and 400 ms, respectively, the third one with UT3-UT4 (62.5 m, 
800 ms), and the last one with UT2~UT4 (89.4 m, 800 ms). Middle: 
Calibrated visibilities of NGC 3783 and geometric ring-model fit (red). 
We derive a ring-fit torus radius of 0.74 + 0.23 mas or 0.16 ± 0.05 pc. 
Bottom: Closure phases plotted versus wavelength. 



errors do not allow us to detect any asymmetry. We also tried 
to derive calibrated visibilities from the 2011 data, but without 
success since the transfer function was unstable. 



To interpret our /T-band visibilities (Fig. [T] middle), we first fit- 
ted a geometric thin-ring model (i.e., ring width - outer radius 
minus inner radius = 0) to the visibilities and derived a ring-fit 
radius of ~0.67 mas. This is only a very rough estimate of the 
torus size since the observed visibilities may not only depend 
on the torus but also on the underlying galaxy within the 60 mas 
field-of-view (FOV) of AMBER and on t he accretion disk (AD), 
which is thought to remain unresolved (Kishimoto et al.' 20071 
|2009a b). Therefore, we have to estimate the flux contributions 
from the host galaxy and the AD point source and take these 
contributions into account when fitting the visibilities. 

From the K-hand image of NGC 3783 in the 2MASS cata- 
log, we estimated the flux contribution of the host galaxy within 
the 60 mas AMBER FOV to 0.5 + 0.2 %. To obtain the NIR flux 
from the torus and the AD, we used two-dimensional fits to sep- 
arate the point-like core component in the 2MASS J-, H-, and 
/T-band images from the underlying host gala xy (see Table |2), 
followi ng the same procedure as described by iKishimoto et al.l 
(I2009ah . Using the derived NIR core fluxes, we can estimate the 
flux contribution of the AD component in the K band. We as- 
sume here that the core component flux originates from the hot 
dust and from the AD. Therefore, we fitted a power-law spec- 
trum for t he AD and a blackbody for the dust emission, as de- 
scribed in IKishimoto et al.l (l2009al) . We also applied a small cor- 
rection for Galactic reddening with Eg-v = 0.119. By assigning 
an uncertainty of the NIR AD spectral index of 0.3, we also ob- 
tained the uncertainty of the ^T-band AD flux contribution. The 
AD flux fraction in the ^T-band was estimated to be as small as 
21 ± 7%, which is similar to the values of several other AGNs 
reported by Kishimoto et al. (2007, 2009a). If we now take into 
account these estimated flux contributions of ~21% from the un- 
resolved AD and of ~Q.5% from the host galaxy, we can derive 
the visibilities of the torus itself and can fit the radius of the torus 
(i.e., this radius is the only fit parameter; the AD contributes just 
a constant of 0.21 to the total visibilities). We derive a torus ra- 
dius /?torus of 0.74 + 0.23 mas or 0.16 + 0.05 pc (thin-ring fit; see 
Fig. [U middle, red curve). 
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Table 2. NGC 3783 torus radius Rtoms, 2MASS fluxes of the nuclear core, and flux contributions of the host galaxy and the AD. 



/flux 


//flux 


KMx 


/ 


H 


K 


host 


AD 


^toms 


^torus 




(mjy) 


(mJy) 


(mJy) 


(mag.) 


(mag.) 


(mag.) 


fraction" 


fraction'' 


(mas) 


(pc) 


(pc) 


18.8 


34.2 


61.8 


12.3 


11.2 


10.1 


0.0()5±().002 


0.21 ±0.07 


0.74 ± 0.23 


0.16 ±0.05 


0.071 ±0.025 



"A'- band flux co ntribution in the AMBER FOV. ''AD flux contribution to the point-like core in the 2MASS A^-band image. ''Reverberation radius 
Rr^ dGlassll 19921) . 
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Fig. 2. /T-band toru s radii of NGC 3783 (red dot) and eight other 
AGNs (blue; from iKishimo to et al. 2011a) versus their UV lu- 
minosities. The bl ack symbols and the dotted line are the rever- 
beration radii Rj,, dGlassll 19921; ISuganuma et al.ll2006l) and tiieir 
fit curve, respectively. 



4. Interpretation and discussion 

4.1. NIR interferometric and reverberation radii 

Figure |2] compares the derived ring-fit radius /Jtoms ~ 0.16 pc of 
NGC 3783 (red) with eight interferometric ^T-band radii (blue) 
reported by Kishimoto et al. (2009a, 2011a). These radii are 
plotted against the UV luminosity L, defined as a scaled V-band 
luminosity of 6 vfv(V), with the V flu x extrapolated from the 
flux at 1.2 fjm (IKishimoto et al.l [20071) . We can compare these 
torus radii with reverberation radii Rj^ (black) derived from the 
light traveling distances coiTes ponding to the time lag between 
the /T-band and the UV/optical (ISuganuma et al.l2006l) . They are 
known to be proportional to ~L'^^ and are likely probing the dust 
sublimation radius. The dotted line is the fit curve of the rever- 
beration radii (different luminosity values are obtained for the 
same object because of variability and uncertainties of the lu- 
minosity derivation). The reverberation radius of NGC 3783 is 
~0.071 pc, which is smaller than the interferometric torus ra- 
dius /?toms ~ 0.16 pc (Sect. 3). Figure |2] shows that several in- 
terferometric torus radii are larger than the reverberatio n radii. 
Our inteipretation is that the inteiferometric torus radii are aver- 
ages over the radial dust distribution that emits the /T-band light, 
whereas the reverberation radii probably trace the dust closer to 
the inner dust torus boundary radius (Kishimoto et al. 2009a). 
Furthermore, we note that Rtoms ~ 0.16 pc is a fit radius calcu- 
lated with a thin-ring model (i.e., ring width = outer radius minus 
inner radius - 0). If a dust distribution with a ceitain thin-iing fit 
radius is ring-like and has a ring width larger than zero, then the 
inner ring radius would be smaller than the thin-ring fit radius. 
We have not fitted a ring model with a larger ring width, since the 
ring width cannot be constrained with the available visibilities. 
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Fig. 3. Temperature/density-gradient model including an addi- 
tional inner hot 1400 K ring component. Top: SED observa- 
tions (black and gray symbols), model SEDs (black solid line: 
model including the 1400 K ring component; black dotted line: 
model without the 1400 K component), and correlated fluxes 
(yellow, green, and blue; see Paper I for more details). The dif- 
ferent colors (see top color bar) correspond to diff'erent spatial 
wavelengths measured in units of the dust sublimation radius Ri^. 
Bottom: New NIR visibilities (purple symbols) and our MIR vis- 
ibilities from Paper I (the symbols with colors from green to red 
coiTespond to 8.5 to 13 /urn; see color-coding bar; note that the 
spatial frequency is given in units of cycles per Ri^). The solid 
and dashed lines are the visibilities of the temperature/density- 
gradient model including an additional inner 1400 K ring. The 
red, green, and purple lines are the model visibilities at 13, 8.5, 
and 2.2 fim, respectively (solid lines: model curves for the PA 
alo ng the eq uatorial axis; dashed line: along the polar direction; 
see lPaper ih . The dotted lines are the visibilities and SED of the 
same temperature/density-gradient model, but without an inner 
hot 1400 K ring component. 
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4.2. Simultaneous modeling of the NIR AMBER visibilities, 
the MIR MIDI visibilities, and the SED 

Mid-infrared (MIR) MIDI inteiferometry of NGC 3783 was re- 
ported by Beck ert et al. (2008), Kishimoto et al. (2009b), and 
iKishimoto et alJ (l201 Ibi) (= iPaper D . For the interpretation 
of these observations, we used a temperature/density-gradient 
model including an additional hot inner ring component with 
a temperature of 1400 K (Paper I). This simple model assumes 
that the face-on surface brightness distribution of the torus is 
dominated by the IR radiation from dust clouds directly illumi- 
nated and heated by the AD. These dust clouds are probably 
located near the torus surface since clouds deep inside the torus 
are not directly illuminated. 

The surface brightness distribution of this 
temperature/density-gradient model depends on two distri- 
butions, namely a radial temperature and a radial surface 
density distribution (see Eq. 8 in Paper I). The maximum dust 
temperature T,^a^{r) at distance r is assumed to be proportional 
to (r/Rinf, where r is the radial distance, j3 is the power-law 
index, and Ri^ is the dust sublimati on radius em pirically given 
by the NIR reverberation radius R^^ (lGlasslll992l) . i.e., we define 

^R^^ (Paper I, Eq. 1). 

Furthermore, the surface brightness distribution depends on 
the surface density function /e(r) = foir/Rin)" of the heated dust 
clouds near the surface (power law with index a). The emissivity 
factor /() is equal to f^(r) at r equal to the sublimation radius 
Our IR observations are only sensitive to the dust clouds near the 
surface, which have the temperature T^.^^(r), and not to the cold 
dust inside the torus. f^{r) can be regarded as a surface filling 
factor multiplied by the emissivity (Paper I). If the emissivity of 
optically thick illuminated clouds does not depend sensitively on 
the radial distance from the illu minating source or the obse rving 
wavelength (e.g., see Fig. 3 of iHonig & Kishimoto! |20 1 Ol) . the 
factor /(r(r) is roughly proportional to the radial surface density 
distribution of the heated dust. 

Interestingly, the application of this temperature/density- 
gradient model to several AGNs in lPaper ll and to the NGC 3783 
observations reported in this paper (see Fig. |3) has shown that 
an additional inner hot model component is required with a tem- 
perature of 1400 K and a radius of one or a few dust sublimation 
radii in order to explain all observations. This hot component 
might play a similar role as the puffed-up inner rim discovered 
in several young stellar objects near the dust sublimation radius. 

We simultaneously fitted this temperature/density-gradient 
model including an inner 1400 K ring t o our new /T-band data as 
well as the MIR data and the SED from lPaper It The goal of this 
modeling is to further constrain physical parameters of the dust 
distribution. Figure|3]shows that this model is able to simultane- 
ously reproduce all NIR and MIR visibilities as well as the SED. 
In Fig.[3](bottom), the NIR visibilities (purple) and the MIR vis- 
ibilities are shown (from green to red, wavelengths 8.5 to 13 fim; 
see color-coding bar). The purple, red, and green curves are the 
model visibilities at 2.2, 13, and 8.5 yum, respectively. The solid 
lines are the model curves along the PA of the equatorial axis, 
the dashed line along the polar direction, as defined by optical 
polarization measurements (see Paper I, where this elliptical and 
a circular symmetric model are presented). 

If there is no hot inner ring component added to the above 
temperature/density-gradient model, then the /T-band model vis- 
ibilities (blue dotted line in Fig. [3] bottom) are systematically 
higher than observed and the model SED (black dotted line in 
Fig. [3] top) has a deficiency in the NIR. Therefore, the above 
inner 1400 K ring component is required to explain the data. 



This new modeling including A'-band visibilities (Fig. O is 
more detailed than that in Paper I. Some of the parameters are 
similar as in Paper I; a temperature power-law index /3 - -0.37, 
density index a = 0.07, and emissivity factor /o - 0.12 for the 
power-law compo nent of t he elliptical model with an inner tem- 
perature of 700 K dPaper J, Table 8). However, in our new model- 
ing, the emissivity factor of the inner 1400 K ring is 0.038+0.016 
and the radius of the hot 1400 K ring is 2.29 + 0.47/?^^ or 
~0.16pc (with theabove /?Tj. = 0.071 pc), which is no longer 
fixed to 1 Rj^ as in lPaper j . This 1400 K ring radius of ~ 2.29/?^^, 
which is relatively large compared to the reverberation radius 
Rt^, is a representative thin-ring radius and not the inner radius 
of an extended ring (see discussion in Sect. 4.1). This large ra- 
dius probably indicates a relatively shallow, extended innermost 
dusty structure (iPaper lb in NGC 3783. 

5. Conclusion 

We have derived a torus radius of 0.74 + 0.23 mas or 0.16 + 
0.05 pc (thin-ring fit). To derive this NGC 3783 torus ra- 
dius, we took into account an estimated relative flux contribu- 
tion of 0.5% from the underlying host galaxy in the 60 mas 
AMBER FOV and 21% from the unresolved accretion disk. 
This torus radius is approximately 2.3 times larger than the K- 
band reverberation radius Rj^, ~ 0.071 pc (see discussion in 
Sect. 4.1). For the interpretation of the observations, we em- 
ployed a temperature/density-gradient model including a hot in- 
ner 1400 K ring. We simultaneously fitted our new N IR AMBER 
visibilities, the MIR MIDI visibilities from'Paper J, and the SED 
to constrain physical parameters of the dust distribution. For the 
power-law component of the model, we derived a temperature 
power-law index /3 ~ -0.37 and a surface density index a ~ 
0.07. For the required 1400 K ring component, a radius of ~2.3 
reverbe ration ra dii or ~0.16pc was found, whereas in the mod- 
eling in iPaper fl. the 1400 K ring radius was assumed to be one 
reverberation radius. This 1400 K ring radius of ~0.16pc, which 
is relatively large compared to the reverberation radius, is a rep- 
resentative thin-ring radius and not the inner radius of an ex- 
tended ring. This radius probably indicates a relatively shallow, 
extended inner du sty stiTicture. Our study of NGC 3783 and the 
results in Paper f show that the simultaneous modeling of both 
NIR and MIR interferometric observations is a powerful tool for 
future detailed studies of AGN tori. 

Acknowledgements. We thank the ESO VLTI team on Paranal for the excellent 
collaboration and the referee for his valuable comments. This pubhcation makes 
use of the SIMBAD database operated at CDS, Strasbourg, France. Some of 
the data presented here were reduced using the publicly available data-reduction 
software package amdlib kindly provided by the Jean-Maiie Mariotti Center 
I http://www.jmmc. fr/data_processing_amber.htm.i. 



References 

Antonucci, R. 1993, ARA&A, 31, 473 

Beckert, T, Driebe, T, Honig, S. F., & Weigelt, G. 2008, A&A, 486, L17 

Burtscher, L., Jaffe, W., Raban, D., et al. 2009, ApJ, 705, L53 

Buitscher, L., Meisenheimer, K., Jaffe, W., Tristram, K. R. W., & Rottgering, 

H. J. A. 2010, PASA, 27, 490 
Chelli, A., Utrera, O. H., & Duvert, G. 2009, A&A, 502, 705 
Glass, I. S. 1992, MNRAS, 256, 23? 
Honig, S. F & Kishimoto, M. 2010, A&A, 523, A27 

Jaffe, W., Meisenheimer, K., Rottgering, H. J. A., et al. 2004, Nature, 429, 47 
Kishimoto, M., Honig, S. F, Antonucci, R., et al. 2011a, A&A, 527, A121 + 
Kishimoto, M., Honig, S. F, Antonucci, R., et al. 2009a, A&A, 507, L57 
Kishimoto, M., Honig, S. F, Antonucci, R., et al. 2011b, A&A, 536, A78, = 
Paper I 

Kishimoto, M., Honig, S. F, Beckert, T., & Weigelt, G. 2007, A&A, 476, 713 



4 



G. Weigelt et al.: VLTI/AMBER observations of the Seyfert AGN NGC 3783 



Kishimoto, M., Honig, S. R, Tristram, K. R. W., & Weigelt, G. 2009b, A&A, 

493, L57 

Kreplin, A., Kraus, S., Hofmann, K.-H., et al. 2012, A&A, 537, A103 

Krolik, J. H. & Begelman, M. C. 1988, ApJ, 329, 702 

Meisenheimer, K., Tristram, K. R. W., Jaffe, W., et al. 2007, A&A, 471, 453 

Petrov, R. G., Malbet, R, Weigelt, G., et al. 2007, A&A, 464, 1 

Pott, J.-U., Malkan, M. A., Rlitzur, M., et al. 2010, ApJ, 715, 736 

Raban, D., Jaffe, W., Rottgering, H., Meisenheimer, K., & Tristram, K. R. W. 

2009, MNRAS, 394, 1325 
Suganuma, M., Yoshii, Y., Kobayashi, ¥., et al. 2006, ApJ, 639, 46 
Swain, M., Vasisiit, G., Akeson, R., et al. 2003, ApJ, 596, L163 
Tatulli, E., Millour, R, Chelli, A., et al. 2007, A&A, 464, 29 
Tristram, K. R. W., Meisenheimer, K., Jaffe, W., et al. 2007, A&A, 474, 837 
Tristram, K. R. W., Raban, D., Meisenheimer, K., et al. 2009, A&A, 502, 67 
Tristram, K. R. W. & Schartmann, M. 201 1, A&A, 531, A99+ 
Wittkowski, M., KerveUa, R, Arsenault, R., et al. 2004, A&A, 418, L39 



5 



